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The astrophysics community is considering plans for a variety of gamma-ray telescopes (including 
ACT, GRIPS, and AdEPT) in the energy range 1-100 MeV, which can fill in the so-called “MeV 
gap” in current sensitivity. We investigate the utility of such detectors for the study of low-mass 
dark matter annihilation or decay. For annihilating (decaying) dark matter with a mass below about 
140 MeV (280 MeV) and couplings to first generation quarks, the final states will be dominated by 
photons or neutral pions, producing striking signals in gamma-ray telescopes. We determine the 
sensitivity of future detectors to the kinematically allowed final states. In particular, we find that 
planned detectors can improve on current sensitivity to this class of models by up to a few orders 
of magnitude. 


I. INTRODUCTION 

The evidence for nonbaryonic, cold dark matter is overwhelming, but there has yet to be any clear indication 
that dark matter interacts with the Standard Model (SM). If such interactions exist, it is possible to obtain indirect 
evidence for dark matter by detecting SM particles that are the products of dark matter annihilations or decays. 
The primary difficulty is discriminating these particles from the astrophysical foregrounds, which themselves are not 
always well understood. Focusing on sharp or distinct spectral features is a straightforward way to enhance the 
signal-to-background ratio and boost detector sensitivity. In particular, a monochromatic gamma-ray line would be a 
reliable indirect detection signal, since it would be difficult to explain such a signal with conventional astrophysics [1]. 

Searching for gamma-ray lines as dark matter signatures is a well-studied topic (e.g., see Ref. [2]). The most re¬ 
cent searches with gamma-ray satellites^ have been performed by the Fermi-LAT Collaboration [3-5]. The upcoming 
space-based telescope GAMMA-400 [6] is set to launch in 2018 and will have an energy range that overlaps with 
Fermi and extends up to 3 TeV. There is, however, a significant “MeV gap” [7] in gamma-ray detector sensitivity 
in the 0.1-100 MeV range. There have been proposals to address this gap with future experiments, such as the Ad¬ 
vanced Compton Telescope (ACT) [8], the Advanced Pair Telescope (APT) [9], the Gamma-Ray Imaging, Polarimetry 
and Spectroscopy (GRIPS) detector [7], the Advanced Energetic Pair Telescope (AdEPT) [10], the Pair-Production 
Gamma-Ray Unit (PANGU) [11], and the Compton Spectrometer and Imager (COSI) [12]. Also, the ASTROGAM 
Collaboration^ is currently working on the research and development of a space-based mission for an MeV instrument. 
Given these efforts by the astrophysics community, it is especially interesting to consider models of low-mass dark 
matter, which can produce distinctive photon signatures in this energy range, particularly since the latest direct- 
detection experiments (such as LUX [13], Super-CDMS [14], and GDMSlite [15]) lose sensitivity for mx GeV, 
making indirect searches all the more important. 

There is a class of dark matter models that will generally produce distinctive photon signatures in the 1-100 MeV 
range: low-mass (-^/s ^ 280 MeV) annihilating/decaying dark matter that couples significantly to first-generation 
quarks. Eor such models, the annihilation/decay products are tightly constrained by kinematics; the only available 
nonleptonic two-body hnal states are 77, 77r°, and These channels will dominate in most of the parameter 

space, and all of these channels produce sharp features, such as lines or boxes, in the photon spectrum. Our focus 
will be on the constraints one may place on these channels from current data and from future telescopes. 

While the annihilation/decay of dark matter to 77 is well studied, the final states 77r° and 7r°7r° are less so. Similar 
studies have investigated monoenergetic photons produced in the processes XX —> ')<fi and X —)■ jcf), with h or Z 
playing the role of (j) [1, 16-18]. Additionally, the cascade process XX cfcj) ^ 7777 produces a box-shaped spectrum 
from (j) decay, and Eermi-LAT data can constrain such channels for dark matter masses above 5 GeV [19]. In our 
case, however, we consider tt*^ in the role of 4>, at a much different energy range. There have also been models of dark 
matter decay that give monoenergetic photons in our energy range of interest [20-22]. 

The plan of this paper is as follows. In Sec. II we introduce the class of low-mass dark matter models that give rise 
to sharp gamma-ray spectral features, and in Sec. Ill we briefly discuss the cosmological bounds on these models. We 


^ Ground-based Imaging Atmospheric Cherenkov Telescopes are also able to perform indirect-detection searches in dwarf spheroidal 
galaxies, but their energy thresholds are ~ 100 GeV, above our energy range of interest. 

^ http://astrogam.laps.inaf.it/index.html. 
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then go into details in Sec. IV on the astrophysical signatures for the indirect detection of dark matter. The specific 
systems we focus on are the gamma-ray diffuse background in Sec. V and dwarf spheroidal galaxies in Sec. VI. We 
comment on possible collider constraints in Sec. VII and conclude in Sec. VIII. 


II. DARK MATTER INTERACTIONS 

We consider low-mass dark matter to be a particle X that couples directly to first-generation SM quarks only. 
Generally, we can consider either dark matter annihilation or decay to SM final states. Since the dark matter initial 
state is assumed to be neutral under all unbroken gauge symmetries of the SM, the only remaining relevant symmetries 
to consider are C, P, and T. The only possible two-body final states that are kinematically accessible are as follows. 

(i) 77 : Accessible at all energies. The final state is C-even. 

(ii) 77r°: Accessible for y's > mj^o. The final state is C-odd. 

(hi) 7r°7r°: Accessible for y/s > 2m^o. The final state is C-even. 

(iv) 7r+7r“: Accessible for y/s > 2m^±. The final state is C-even or C-odd. 

(v) U {£ = e,fi, v): Accessible for y/s > 2m(. The final state is either C-odd or is weak suppressed. 

Note that if X decays rather than annihilates to SM particles, then these final states are only allowed if A is a boson. 
If we assume that weak interactions, which are suppressed by a factor sGp, are negligible, then the only leptonic states 
which can be produced are the C-odd e+e” and final states via an intermediate off-shell photon that couples to 

a quark loop. But the annihilation/decay rate in this channel will be suppressed by a factor so it will be dominant 
only if no other channel can be produced. Since the photon spectrum arising from this channel is not as distinctive as 
the others, we will ignore it from here on. Three-body final states are also accessible if additional photons or neutrinos 
are emitted; however, these processes come with additional suppression factors of a or sGf, respectively, and can 
thus also be ignored unless all of the channels above are heavily suppressed. It is worth mentioning that final-state 
radiation (FSR) and virtual internal bremsstrahlung (VIB) are capable of producing somewhat distinctive features in 
the photon spectrum. FSR yields a continuous XjE spectrum [23], distinct from the astrophysical foregrounds that 
have a spectral index around —2 (see Sec. VB). It may be possible to observe the upper edge of such a spectrum, 
determined by mx for annihilation or mxl2. for decay [24]. Similarly, VIB can produce a broad bump near ^ mx 
(for annihilation), which mimics a line in a detector with poor energy resolution, but the details of the spectrum tend 
to be more model dependent [25]. 

For y/s > 2mTr± it is kinematically possible to produce the 7 r“'' 7 r“ final state. As with charged leptons, FSR can 
produce a sharp edge in the photon spectrum; however, this process could potentially compete with the 77 r° or 
TT^TT^ channel, diminishing the prospect of observation. Moreover, even though these states with one or two neutral 
pions could be produced with non-negligible branching fractions, their photon spectrum becomes very broad if y/s is 
significantly greater than 2 m,r 0 and thus less interesting from the point of view of detectability at future experiments. 
Additionally, one would expect the branching fraction to the 77 final state to be suppressed by ~ c? relative to the 
7 r+ 7 r“ state. We thus see that this class of models produces the most interesting photon spectrum, from the point of 
view of detectability, for the range y/s ^ 2mT^±. We limit our study to this mass range. 

If we assume that weak interactions are negligible, then the only potential source of C violation is via the coupling 
of dark matter to the first-generation quarks. If we insist that this coupling is instead C invariant, then one may 
classify allowed final states by the C quantum number of the initial state. For example, in a C-invariant effective 
field theory, a C-odd initial dark matter state can only produce the 771 ° or U two-body final states; if the 777 ° state 
is kinematically allowed, then it will dominate. An example of this type of fundamental interaction (if A is a Dirac 
fermion) would be the interaction (l/M^)(A 7 ^A)(g 7 ^g), where q = u,d. In this case, the AA initial state with a 
nonvanishing matrix element is C odd. In the low-energy effective field theory involving only A, 7 , and 7 r° (excluding 
leptons), the lowest-dimension Lagrangian interaction one can write that is Lorentz and C invariant and linear in 
A 7 '"A is 




^(A7^A)F'^'’(5"7rO)e^,.,,, . 


(1) 
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The coefficient A can be related to the one-loop diagram mediating the tt^ —t 77 decay interaction (by replacing one 
photon with the dark matter vector current), which is determined by the chiral anomaly.^ We would then expect 
1 /A^ ^ (e/167r^)(l/M^/,r)j where is the pion decay constant. 

In general, however, the dark matter-SM quark-level interaction need not preserve C. As such, dark matter 
annihilations (decays) are able to produce all three final states of interest—provided they are kinematically accessible— 
with branching ratios Ajr, and A-^ 'D-ym and V-y), respectively. These branching ratios depend on the specific 
UV model, but we leave them here as parameters to keep our analysis general. 


A. Photon spectra 


The prompt photons simply have i5-function spectra in the dark matter center-of-mass frame: 


dN^ 

~dE 


5{E - Eo) , 


( 2 ) 


for each photon produced with an energy Eg- Any tt^ subsequently decays to two secondary photons with a branching 
ratio of ^ 99% [26] . The decay is isotropic in the 7 r° rest frame, and boosting to the dark matter center-of-mass frame 
results in a box-shaped photon spectrum [27] 


ItE 


2 
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[Q{E-E_)-Q{E-E+)] , 


( 3 ) 


where E± are the kinematic edges and /S.E = E+ — E- is the box width. Thus, the annihilation/decay processes 
yielding prompt photons or neutral pions produce gamma spectra with sharp features. We summarize the kinematics 
below. 


(i) 7 r*^ 7 r*^: The photon spectrum is box shaped, given by twice that in Eq. (3), with kinematic edges and box width 


E± = 





( 4 ) 


(ii) 77 r°: The prompt photon produces a line distribution, given by Eq. (2), with energy 


E„=^ll- 


The spectrum from the pion decay is given by Eq. (3), with kinematic edges and box width 


E+ = 
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(iii) 77 : The photon spectrum is a line, given by twice that in (2), at the energy 


( 5 ) 

( 6 ) 

( 7 ) 


If the box spectrum is very narrow, a detector will not be able to resolve the box shape and will instead observe 
a signal that is indistinguishable from a line. Since the box width is larger for a more highly boosted pion, a pion 
produced nearly at rest ( 1 / 5/2 « 771,^0 for the 7r°7r° channel and ^/s « 777,^0 for the 777 ° channel) produces two photons 
with energies close to mjroj2. At the upper end of the kinematic range we consider (y^ = 2777,r±)i tho width of the 
box spectrum for the 777 ° and channels is ^ 106.9 MeV and ~ 35.5 MeV, respectively. 


® The Xjifi A coupling is not related via isospin to any potential coupling of dark matter to 7r^, because electromagnetic interactions 
violate isospin near-maximally. 














4 


III. COSMOLOGICAL BOUNDS 

Searches for particle dark matter are typically performed by looking for signals from production at colliders, direct 
detection, or indirect detection. Collider constraints are highly model dependent, and later in Sec. VII we will briefly 
discuss bounds in the context of the example model presented in Sec. II. But since we are interested in the broader 
class of models, we do not consider collider constraints in our overall analysis. Current direct-detection experiments 
are sensitive to dark matter masses above ~ 1 GeV, so their results do not apply to dark matter at much lower 
masses. However, the annihilation/decay of light dark matter can be constrained by its effect on the photon spectrum 
observed by indirect-detection experiments, which is the main focus this paper. 

Another constraint arises from the effects of dark matter annihilation/decay on cosmological history. Before in¬ 
vestigating the indirect detection of dark matter, we first discuss these limits from the early Universe. Dark matter 
annihilations and decays can inject energy into the primordial plasma, altering the predictions of big bang nucleosyn¬ 
thesis (BBN) [28-33] and the spectrum of the cosmic microwave background (CMB) [34-45]. For dark matter masses 
0(100 MeV), CMB limits for s-wave annihilation are more stringent than those from BBN [33]. If we consider the 
latest Planck limits [46], the upper bound on the thermalized annihilation cross section is 

{av) < 3 X 10-28 cmVs • (8) 

Here we assume that dark matter annihilation produces only photons, since the neutral pions that are produced rapidly 
decay into photons on a time scale less than the Hubble time scale during recombination. This limit is about an order 
of magnitude stronger than the previous WMAP9 limit [42]. If the annihilation process is instead p-wave, the CMB 
limits weaken [45] and present-day limits from the diffuse gamma-ray background can dominate [47]. In this case, 
future gamma-ray telescopes providing better measurements of the diffuse background can also better constrain dark 
matter p-wave annihilations. The more challenging question is whether or not future telescopes will have sensitivities 
to reach the strong bound in Eq. (8) , so here we focus on the simpler analysis of s-wave annihilation. For decaying dark 
matter, we use the numerical tools presented in Ref. [41] to analyze the case of dark matter decaying to monoenergetic 
photons, and the resulting lower bound on the lifetime from Planck is r > s. The annihilation and decay 

bounds are displayed as black lines in the left and right panels, respectively, of Fig. 1. 

Note that the most stringent current constraints on dark matter annihilation come from its effect on the CMB, 
while CMB constraints on dark matter decay are less stringent than those arising from the observed diffuse gamma-ray 
background, as we discuss in Sec. V A. CMB constraints are much more severe for the case of dark matter annihilation, 
because the annihilation rate is proportional to the dark matter density squared, while the decay rate is proportional 
to just a single power of the density; the annihilation rate thus has a large enhancement in the early Universe relative 
to the current epoch. 


IV. GAMMA-RAY FLUX 


The differential flux of gamma rays from annihilating dark matter is 


^ J (av) 

dE dfl 47r 2 to^ * 


dE' 


Re{E 


E') dE' , 


( 9 ) 


where mx is the dark matter mass, {av) is the velocity-averaged annihilation cross section, and is the smearing 
function due to the limited detector energy resolution ue = iE. Here we assume that dark matter is its own 
antiparticle; if it is not, then the expression for the flux is rescaled by a factor of 1/2. We estimate the detector 
smearing function with a Gaussian, which must be convolved with the true spectrum. For a gamma ray with true 
energy E', the probability that it will be reconstructed with an energy E is [48, 49] 


Re{E-E') 


1 

-x= -exp 

V^eE' 


{E-E')^\ 
2e‘^E'^ ) 


( 10 ) 


As a result, an injected monochromatic line distribution will appear at the detector as a Gaussian distribution, 
centered at the line energy Uq, with a half-width given by eEo. The kinematic edges of the pion decay box-shaped 
spectrum will be smeared out as well. The sum in Eq. (9) runs over all dark matter annihilation processes i with 
branching ratios Ai and spectra dN^/dE. The total gamma-ray spectrum has contributions from continuum spectra 
and from monochromatic lines, and we are only concerned with spectra that are isotropic. 
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Detector Source Energy Range [MeV] e PSF Agfi [cm^] 


ACT 

[8] 

0.2-10 

1% 

1° 

1000 

GRIPS 

[7] 

0.2-80 

3% 

1.5° 

200 

AdEPT 

[10] 

5-200 

15% 

0.5° 

600 

COMPTEL 

[59, 60] 

0.8-30 

2% 

2° 

50 

EGRET 

[61] 

30 MeV-10 GeV 

12.5% 

2.8° 

1000 

Fermi-LAT 

[62] 

20 MeV-300 GeV 

7.5% 

2° 

4000 

GAMMA-400 

[6] 

100 MeV-3 TeV 

12% 

2° 

3000 


TABLE I: Experimental parameters used to determine indirect-detection bounds for dark matter annihilation and decay. The 
PSF and AbB values are not needed in the analysis for COMPTEL, EGRET, and Fermi, but they are included for comparison. 

All of the terms in Eq. (9) depend only on the particle properties of the dark matter, except for the astrophysical 
J-factor 

JW = [ dl , (11) 

Jlos 

where the integral is taken over the line of sight, and p(r) is the dark matter energy density distribution. The 
coordinate r is centered on the dark matter halo of interest, which is a distance D from the sun. The quantity '0 
is the angle between the line of sight and the line connecting the sun with the center of the halo. Thus, we have 
r — {D^ — 2Dl cosip + , and the line-of-sight integral is from 0 to ^max = — siiP ipD^yP + D cosip [50]. To 

find the total dark matter flux from a particular region of the sky, we need the J-factor integrated over that region: 

j = Jdn^j{iP). ( 12 ) 

The total number of photons with energies between Ei and E 2 detected from a source with an integrated J-factor 
given by observed over a period of time Tobs, is 

= T,bsAeS<^ = ^obsJteff^ ^ ^ 

where Aes is the effective area of the detector. The analogous expressions for decaying dark matter are obtained by 
making the substitutions {av) /2mx —> E and Ai —>■ Vi in Eqs. (9) and (13), and —>■ p in Eq. (11), where T is the 
decay width and the sum runs over decay processes i with branching ratios Vi. 

There are a few ideal systems to search for an indirect-detection signal from dark matter: dwarf spheroidal galaxies, 
the diffuse gamma-ray background, and the Galactic center (GC). We do not focus on regions near the GC in this 
work, due to the large astrophysical uncertainties involved. For example, the uncertainty in the spectrum arising from 
astrophysical backgrounds is very difficult to quantify. In addition, the central halo profile is uncertain due to effects 
from a supermassive black hole [51-54], from interactions with baryons [55, 56], and from baryonic infall [57, 58]. 

In the following sections, we describe the properties of the other two systems and their potential for producing 
signals from dark matter. We treat the three dark matter annihilation/decay channels of interest separately, assuming 
in each case that the branching ratio is 100%. For each of these channels, we determine limits from the diffuse 
background for previous and existing experiments (COMPTEL, EGRET, and Fermi), and we give projected limits 
from both the diffuse background and dwarfs for future experiments (ACT, GRIPS, AdEPT, and GAMMA-400). 
The general detector properties we use are summarized in Table 1. We use more conservative estimates of energy and 
angular resolution (in the most relevant energy regions) provided by the experiments. Actual detectors have complex 
response functions (e.g., ge is a more complicated function of energy, and Aes depends on energy and detected 
angle of incidence) and are typically optimized in the central region of their energy window. Given that we are 
mainly interested in the potential of future experiments rather than setting the most precise constraints on existing 
experiments, we use these baseline numbers for a more straightforward analysis. 


V. INDIRECT DETECTION: DIFFUSE GAMMA-RAY BACKGROUND 

One way to search for dark matter is observing the diffuse gamma-ray background. There are two components 
to the diffuse background: Galactic and extragalactic. We expect the diffuse background to comprise signals from 
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both dark matter and astrophysical processes, but the astrophysical contribution cannot be estimated entirely from 
data: one must instead make some assumptions about an underlying astrophysics model. Additionally, there will be 
enhancements for Galactic and extragalactic dark matter annihilation due to substructure [63-65]; however, different 
models can yield a wide range of results [4], so we ignore substructure enhancements and obtain more conservative 
bounds on dark matter annihilation. 

The Galactic contribution from dark matter is simply given by Eq. (9). We use publicly provided calculations of 
[ 66 ] with an Nevarro-Frenk-White profile; far from the GG, the J-factor becomes less sensitive to the details of 
the central profile, and the integrated J-factors for various profiles differ by 0(1) factors at high Galactic latitudes [ 66 ]. 
For the extragalactic background, we assume the dark matter contribution is isotropic. Photons detected with energy 
Eq were emitted at redshift z with energy E{z) = (1 + z)Eq. Although redshifting has the effect of smearing out an 
otherwise sharp signal from a monoenergetic photon or from neutral pion decay, the extragalactic flux is subdominant 
to the Galactic flux, particularly in the case of annihilation. Photons emitted at redshift z 1 are relevant, because 
they provide a small enhancement on top of the Galactic signal, within the energy resolution of the detector. 


A. Limits from existing data 

For the energy range of interest, COMPTEL, EGRET, and Fermi-LAT give the most relevant data of the differential 
flux for the Galactic and extragalactic spectrum'* at high Galactic latitudes b. COMPTEL [59] observed the region 
15] > 30°, EGRET [67] observed 20° < [fe] < 60° (see their Fig. 8 , region E), and Fermi [ 68 ] observed [fe] > 20° (see 
their Fig. 4). Note that although Fermi is capable of detecting photons down to energies of 20 MeV, the data for 
the diffuse analysis only goes down to energies of 100 MeV in order to sufficiently reduce the cosmic-ray background 
produced in the Earth’s atmosphere. 

To determine a conservative limit on the lifetime (cross section) for decaying (annihilating) dark matter, we follow 
the procedure from Ref. [47]: we require that the number of events expected from a dark matter signal in each energy 
bin does not exceed the observed number of counts by 2a. The resulting limits are shown in Fig. 1 for each of the 
three channels of interest: 77 (solid), 771 ° (dashed), and (dotted). Prompt photons will typically fall within one 
or two energy bins (depending on the effects of the energy resolution and the bin placement relative to the photon 
energy). Photons arising from 7 r° decay will fall roughly uniformly within the the energy range bounded by and 
A-|_; the most important energy bin for these photons is the highest energy bin within this range, since the observed 
flux falls as a function of energy. 

The effect of the discrete binning is clear, particularly for the 77 channel, and the rounded shape for each bin is due 
to the smearing effects of the detector energy resolution®. For the 777 ° channel, the projected sensitivity of GOMPTEL 
is dominated by the monoenergetic photon and is best at ^/s ~ 160 MeV. Just above threshold in either the 77 or 
77 r° channels, COMPTEL’s sensitivity drops rapidly because any monoenergetic photon is below COMPTEL’s energy 
range, while any photons arising from pion decay are nearly monoenergetic and are above GOMPTEL’s energy range. 
For either the 777 ° or channel, GOMPTEL’s sensitivity to the pion is very weak, since the box feature in the 
spectrum must be very wide to overlap GOMPTEL’s energy range. For the 77 channel at large ^/s, GOMPTEL’s 
sensitivity is dominated by the extragalactic diffuse component, since the gamma rays at low redshift have energies 
above COMPTEL’s energy range. The limits from EGRET and Fermi are much simpler, since m .,^0 lies within their 
energy ranges. 

It is worth noting that a conservative GC analysis for COMPTEL and EGRET [69] is able to set limits on dark 
matter annihilation to 77 that are comparable to the high Galactic latitude bounds in Fig. 1. There is another 
GC analysis for EGRET [70] that constrains the annihilation cross section to 77 at the level (av) ^ 10“®° cm^ for 
^/s = 100 MeV. However, this analysis presupposes that the photon spectrum arising from astrophysical background 
is smooth and can be determined with subleading systematic uncertainties by fitting to data; the allowed dark matter 
contribution arises only from a 2a downward statistical fluctuation in the background in any bin. This analysis is thus 
not as conservative as either the conservative analysis approach used above (and used by Ref. [69]) or the optimistic 
analysis described in Sec. VB, in which systematic uncertainties dominate. 


^ Although COMPTEL presents results for the extragalactic diffuse spectrum, the given measurements do not have the Galactic spectrum 
or point sources removed due to large systematic uncertainties [59]. 

® The sizes of the bins chosen by each experiment are larger than the minimum sizes determined by the energy resolution. 
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Annihilation Constraints 



^fs [MeV] 

Annihilation Constraints 



Decay Constraints 




FIG. 1: CMB and diffuse constraints on dark matter annihilation (left) and decay (right). The top and bottom sets of panels 
display the same information, plotted on a log and linear x axis, respectively. The strongest CMB bound comes from Planck, 
shown in black for annihilation [46] and decay [41]. Conservative bounds from COMPTEL, EGRET, and Fermi for diffuse 
emission are shown in blue, red, and green, respectively, for the channels 77 (solid), 77 r° (dashed), and (dotted). Optimistic 
bounds, described in Sec. VB, can be obtained by improving the COMPTEL, EGRET, and Fermi bounds by a factor of ~3-5, 
9, and 6 , respectively. Vertical dotted lines show the kinematic thresholds at m^o and 2m^o, and is plotted up to m^±. 


B. Limits from future experiments 

For the existing experiments, we used their energy binning for their data. In analyzing future experiments, we use 
an optimized binning strategy in an effort to obtain the best sensitivity. Practically, doing so avoids the discreteness 
(the small dips) of the curves in Fig. 1. The relevant energy bins are the ones encompassing the monoenergetic photon 
and/or the box feature in the photon spectrum arising from pion decay. The bin for the prompt photon is one centered 
at its true energy Eq and has a width 2eEQ (i.e., the bin detects 68 % of the line signal). For the box spectrum, the 
energy bin that gives the best constraint is the one at the upper edge of the box, so we define a bin with an upper 
edge at E+. In either case, if the signal leaks beyond the lower or upper detector threshold, the bin edge is set to 
be at the threshold. For either the 77 or channel, the only energy bin to consider is the one corresponding to 
the line or upper box edge, respectively. For the 77 r° channel, the more constraining bin (if it is inside the detector 
energy window) is typically that for the line, as we saw in Sec. V A. 

The total integrated flux is over the angular coverage of the system of interest and over the energy bin that best 
constrains the dark matter signal. We consider sensitivities which can be obtained with two different strategies. The 
conservative sensitivity (as in Sec. V A) can be made by assuming that diffuse emission may be due solely to dark 
matter, excluding models that would be expected to produce a number of events in excess {2a) of that measured 
in any bin. An optimistic sensitivity can be derived by assuming that the energy spectrum of the astrophysical 
background exhibits no sharp features and can be determined by fitting to the data. One then excludes models that 
would produce an excess of events in any bin even if the astrophysical background in that bin were overestimated by 
a 2a systematic uncertainty. We assume up to a 15% systematic error in the determination of the smooth component 
of the differential flux, which is quoted for EGRET [67]. Note that we need not consider statistical fluctuations for 
the optimistic analysis, since the systematic uncertainties dominate. 

The background is estimated from a single power-law fit to the COMPTEL data [59] over energies 0.8-30 MeV and 


































to EGRET data [67] over energies 30 MeV-10 GeV: 


= 2.74 X 10 ^ f ^ cm ^sr ^MeV ^ . (14) 

dEdn \MeV J ^ ' 

(Fermi data is not used, since it overlaps and agrees well with EGRET data at energies less than a few hundred 
MeV.) The COMPTEL data corresponds to observations over the region | 6 | > 30°, while the EGRET data was an 
average of two regions (20° < | 6 | < 60° and 60° < | 6 | < 90°), weighted by their relative spans of the sky. Note that 
the data from EGRET encompass a slightly larger region compared to GOMPTEL, so the fit underestimates the flux 
at lower energies and overestimates the flux at higher energies. For E < 150 MeV, the EGRET data in the region 
20° < | 6 | < 60° is only a factor of 1.2-1.4 larger than that in the region 60° < | 6 | < 90°. Thus, we use this fit as a 
rough guide for predicting the diffuse background a future telescope would detect. 

The resulting sensitivities are shown in Fig. 2 for observing the region | 6 | > 30°. For either the conservative or 
optimistic analysis, it is important to note that the sensitivity presented here is not determined by the exposure; the 
maximum differential gamma-ray flux which can be attributed to dark matter is already measured by GOMPTEL, 
EGRET, and Fermi, and we assume statistical uncertainties are subleading. Instead, the sensitivity is determined 
by the experiment’s energy resolution. Gertainly, more reliable astrophysical modeling in the future will be able 
to significantly improve sensitivity estimates and constraints on dark matter. In particular, anisotropies in the 
extragalactic diffuse background produced from dark matter annihilations are predicted to be different from those 
produced by standard astrophysical sources [71-74]. However, in the meantime, we translate our ignorance of the 
functional form of the diffuse background into the 15% systematic error previously discussed. 

Comparing the results in Fig. 2 to those in Fig. 1, projected sensitivities from GRIPS are better than all present 
diffuse bounds for 77 and . For 77 , the monoenergetic photon line lies within the energy range of GRIPS up 
to ~ 160 MeV, and the sensitivities for {av) /mx and the lifetime are approximately constant in ^/s. Similarly, 
for TT^TT^, the box spectrum of the pion is narrow enough near i/s ^ 270 MeV to mimic a line signal, but contributes 
less to the sensitivity as the box widens at higher i/s. For jTr^, the pion box provides the dominant contribution 
to sensitivity until it extends too far above the GRIPS upper threshold, for ^/s > 160 MeV, at which point the 
photon line gives the better constraint. GRIPS’s sensitivity will exceed present bounds until it suffers a sudden loss 
of sensitivity for ^/s ^ 240 MeV, where the photon line is above the GRIPS energy threshold; meanwhile, the photon 
line is still inside the energy range of EGRET and Fermi. 

AdEPT sensitivities are comparable to present diffuse bounds, and it is less sensitive than GRIPS, due to its larger 
energy resolution. The exception is in the 77 r° channel for -y/s ^ 240 MeV, where the photon lies within the AdEPT 
energy window but is above the GRIPS upper energy threshold. In fact, for the mass range plotted, all spectral 
features from dark matter annihilation/decay are contained in the AdEPT energy window for all three channels. For 
77 r°, the pion box provides the better sensitivity for ^/s < 200 MeV, above which the line and box signals overlap. The 
energy bin that contains the line also contains part of the box signal, so this bin becomes the more constraining one. 
For 77 , the sensitivities for (av) jmx and the lifetime are approximately constant in except when i/s 10 MeV, 
for which the photon line drops below the AdEPT energy threshold. 

ACT sensitivities are the most strict in its relevant mass range for the 77 channel, with its conservative sensitivity 
lying on the CMB limit for annihilation. Like GRIPS and AdEPT, the sensitivities for (av) jrax and the lifetime are 
approximately constant in ^/s. For the 77 r° channel, as with COMPTEL (but unlike GRIPS and AdEPT), sensitivity 
is weak near the threshold at ~ 135 MeV, because the pion box is above ACT’s energy range, and the low-energy 
photon line is surrounded by a large background. ACT barely does better than EGRET near ^/s ~ 145 MeV, until 
the photon leaves the ACT energy range at larger ^/s. For the mx range of interest, the box spectrum from 7 r° 7 r° is 
at much higher energies than the upper energy threshold of ACT; thus, no ACT limits are included for this channel. 

The sensitivity of GAMMA-400 is on par with Fermi for the 77 channel. The photon line from 77 r° is never 
above GAMMA-400’s lower energy threshold in most of the plotted mx range, so its sensitivity is mainly due to 
the pion box spectrum. However, the true upper edge of the pion box only surpasses GAMMA-400’s lower energy 
threshold for ^/s > 200 MeV; for smaller -y/s, GAMMA-400’s sensitivity arises from the smearing of the box due to 
the energy resolution. Since Fermi has a better energy resolution than GAMMA-400, as quoted in Table I, there is 
more leakage into the energy range of GAMMA-400 for y/s < 200 MeV, resulting in the greater projected sensitivity 
for GAMMA-400. But this result strongly depends on the form of the detector response function and may not hold 
if the actual response function is significantly different from the one we use [Eq. (10)]. The situation for the 7 r° 7 r° 
channel is similar, except that there is no plotted region for which the true pion box lies within GAMMA-400’s energy 
window—its sensitivity is due solely to the smearing of the spectrum from the detector energy resolution. 
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C. Diffuse background in the 0.3—10 MeV range 

We note that the explanation of the origin of the extragalactic background in the MeV region is incomplete 
and remains under investigation. For energies ^ 0.3 MeV, the diffuse background can be largely explained by the 
contribution of active galactic nuclei and Seyfert galaxies [74-77], while for energies > 10 MeV it is believed that 
blazars [78-80] and radio and star-forming galaxies [81, 82] can provide an adequate explanation for the observed 
background. There is also a small but insufficient contribution from Type la supernovae [83, 84] near 1 MeV. As a 
result, there is no clear explanation for the observed (sharply falling) spectrum in the intermediate range 0.3-10 MeV. 

The class of dark matter models in Sec. II can contribute to the observed photon flux in this energy range. However, 
the spectral signatures of redshifted gamma-ray lines and boxes alone are not suitable to explain the sharply decreasing 
intensity of the observed background. It would be interesting to consider whether these distinctive features from dark 
matter annihilation/decay, in combination with contributions from astrophysical sources, could adequately explain 
the observed photon spectrum in this energy range. A more detailed analysis is beyond the scope of this work. 


VI. INDIRECT DETECTION: DWARF SPHEROIDAL GALAXIES 

The Milky Way satellite dwarf spheroidal galaxies are dark matter dominated and thus very faint, so we do not 
expect any significant source of photons from dwarfs due to standard astrophysical processes, making dwarfs a very 
clean set of systems to study. Although the expected signal is low, the dwarfs lie away from the Galactic plane, so 
the astrophysical backgrounds are much smaller than they would be otherwise. Importantly, one can estimate the 
astrophysical background from data by looking slightly off axis. Table II shows the integrated J-factors for 20 dwarfs, 
from Ref. [85]. The total flux from a dwarf is contained within an angle 0max, set by the outermost observed star 
in the halo, from the center of the dwarf. We consider regions 9 > to give zero contribution to the density 

profile, resulting in a conservative value for the integrated J-factor. In determining the dark matter flux measured by 
a detector, the angular integration of the J-factor extends over the region defined by the halo or by the detector point 
spread function (PSF) 68 % containment radius, whichever of the two is larger. We choose Draco as our case study for 
two main reasons: it has large integrated J-factors for both annihilation and decay, and independent calculations of 
J agree well with one another [85-87]. Comparison with other dwarfs is straightforward using Table II; in particular, 
for a given observed flux, the corresponding dark matter annihilation cross section or decay width scales inversely 
with the J-factor in Eq. (9). Performing a stacked analysis from the observation of many dwarfs will further improve 
sensitivity, though such a detailed study is beyond the scope of this work. 

There are recent constraints on dark matter annihilation from searches in dwarf galaxies using Fermi data. Refer¬ 
ence [ 88 ] performed an analysis for annihilation to 77 , and the limit they derived is comparable to CMB limits (see 
their Fig. 11), but the analysis was for mx > 1 GeV and thus it is not included here. 


A. Limits from future experiments 

For analyzing Draco, the relevant energy bins are all the ones encompassing the dark matter signal, whether they 
are from the monoenergetic photon or from the box spectrum from pion decay, since there are relatively low statistics 
expected from dwarf observations. We again use an optimum bin for the photon line, as discussed in Sec. VB, and 
we define a single large “bin” for the box that spans energies from (1 — e)A_ to (1 -I- e)E+, unless detector energy 
thresholds truncate the box spectrum. To find the total integrated flux, the angular integration is over the PSF of 
the detector. The integrated J-factor for a PSF > 1° is very close to simply using 0max instead. Thus, despite Draco 
having one of the largest extended profile emissions [85], we may still treat it as a point source. 

We assume that the only relevant background is the diffuse emission, and unresolved point sources contribute 
negligibly. Any diffuse contribution from dark matter is already incorporated in the data, so the fit (14) represents 
the total background for observing dwarfs. If the expected number of events for a model in relevant energy bins 
(including signal and background) is /i, then the probability of the model not yielding an excess of events above 
background is 


^obs -I 

P{n < nobs) = V] , (15) 

n\ 

n—Q 

where nobs is the number of observed events, estimated by Eq. (14). We then find the 95% exclusion limits for the 
annihilation cross section or decay rate for five years of running. The limits are plotted as hatched regions in Fig. 2. 
The band thickness represents the Icr systematic uncertainty for the integrated J-factor in Table II. 
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Dwarf 

longitude 1 latitude b ffmax 

loglO ^ann 



[deg] 

[deg] 

[deg] [GeV’^cm-®] 

[GeVcm-2] 

Bootes I 

358.08 

69.62 

0.47 

18.24t°;« 

17.90t°;2® 

Carina 

260.11 

-22.22 

1.26 

17.92t[!;}® 

IS.lbtlfs 

Coma Berenices 

241.89 

83.61 

0.31 

19.02t°jJ 

17.961°;^“ 

Canes Venatici I 

74.31 

79.82 

0.53 

17.44t°fs 

17.57t°of3 

Canes Venatici II 

113.58 

82.70 

0.13 

17.65t°:« 

l&.97+_°olt 

Draco 

86.37 

34.72 

1.30 

19.05«i? 

18.97lg;il 

Fornax 

237.10 

-65.65 

2.61 

17.84t°;“ 

17.99lg;“ 

Hercules 

28.73 

36.87 

0.28 

16.86t°:^^ 

16.66lg;« 

Leo I 

225.99 

49.11 

0.45 

17.84t°;?° 

17.9llg;i® 

Leo II 

220.17 

67.23 

0.23 

17.97tUI 

17.241°;! 

Leo IV 

265.44 

56.51 

0.16 

16.32ti;^ 

16.121?:! 

Leo V 

261.86 

58.54 

0.07 

16.37«:®^ 

15 . 871 ?:! 

Leo T 

214.85 

43.66 

0.08 

17 11 

^ ' •-*--‘--0.39 

16.481°;! 

Sculptor 

287.54 

-83.16 

1.94 

18.57t°:°^ 

18.471?;! 

Segue 1 

220.48 

50.43 

0.35 

19.36t°;f5 

17.991?;?? 

Segue 2 

149.43 

-38.14 

0.19 

iD.Zi_Q gg 

15.891?;?? 

Sextans 

243.50 

42.27 

1.70 

17.92«;®® 

18.561?;?? 

Ursa Major I 

159.43 

54.41 

0.43 

17.87t°®3®3 

17.611?:?? 

Ursa Major II 

152.46 

37.44 

0.53 

19.42t°;« 

18.391?;?? 

Ursa Minor 

104.97 

44.80 

1.37 

18.95t°:f8 

18.131?:?? 


TABLE II: Integrated J-factors for dwarf galaxies from Ref. [85]. The integration is centered on the dwarf and extends out 
to an angle of 0max, the point at which the halo density profile is truncated. The la errors are due to the uncertainty in the 
precise shape of the profile. The location of the dwarfs are given by their Galactic coordinates from the NASA/IPAC 

Extragalactic Database (http://ned.ipac.caltech.edu/). 


As expected, dwarfs are the more advantageous target if dark matter annihilates, since the annihilation rate scales 
as the square of the dark matter density; the large density of dark matter within a dwarf then provides a large 
enhancement to the photon production rate. In fact, ACT, GRIPS, and AdEPT sensitivities of the annihilation 
cross section from Draco are able to reach the bound from Planck’s CMB observations, except when GRIPS loses 
sensitivity near i/s ~ 240 MeV in the 771 ° channel. Additionally, AdEPT has a comparable or better sensitivity than 
GRIPS in all three channels, despite its sensitivity being worse for the diffuse background. This result reflects the 
fact that counting statistics are the determining factors for the dwarf sensitivities, giving dwarfs the advantage over 
the systematics-dominated diffuse background. AdEPT is helped by its larger Aeff and by its larger energy window, 
compensating for its poorer energy resolution. 

On the other hand, the dark matter decay rate only scales as one power of the density, so the sensitivities derived 
from a dwarf versus diffuse analysis should not differ as greatly as they do for annihilation. Indeed, the diffuse and 
dwarf decay sensitivities are comparable for GRIPS, and the dwarf sensitivity is only slightly better for AGT, AdEPT, 
and GAMMA-400. 


VII. A COMPARISON TO COLLIDER-BASED BOUNDS 

The interactions of low-mass dark matter with first-generation quarks can be tightly constrained by “monoanything” 
searches [89-94] at the LHC [95-100]. However, the sensitivity of these searches depends on the details of the dark- 
matter-quark effective interaction; it is thus difficult to compare such constraints to the results of the analysis 
presented here, which only depends on the annihilation/decay branching fraction to each channel. However, one 
can gain some insight into the relative sensitivity of indirect and collider search strategies for a particular model, 
namely, fermionic dark matter which couples to first-generation quarks via an effective contact operator, given by 
O = (l/M^)(XrX)(qj^g). 

If the effective contact operator approximation is valid, then data from ATLAS require M > 800 GeV for low-mass 
dark matter [97]. This coupling permits dark matter annihilation to the 771 *^ final state via the effective operator 
given in Eq. (1). If we require that for ^/s = 2m^o this cross section obey {<jav) 10“^ pb (roughly the optimistic 




11 





10 "^“ —^^^—'—'—^^^^^^^^^^ 

270 272 274 276 278 

/s [MeV] 



50 100 150 200 250 


^fs [MeV] 




FIG. 2: Projected sensitivities on dark matter annihilation (left) and decay (right) for channels 77 (top), 77 r° (middle), and 
(bottom). Planck (black) bounds are shown from Fig. 1. The cyan, blue, red, and green shaded regions show the areas 
between the conservative and optimistic sensitivities for the diffuse background for ACT, GRIPS, AdEPT, and GAMMA-400, 
respectively. The hatched regions show the Icr uncertainty bounds for Draco after five years of observation. Vertical dotted 
lines show the kinematic thresholds at m^o and 2m^o (plotting begins at m^o for the 771 ® channel and at 2m^o for the 
channel), and y/s is plotted up to m^± for all three channels. Note that the vertical axes cover a different range for 77 than 
for the other two channels. 


sensitivity of GAMMA-400 searching for a signal from Draco), then one hnds M ^ 10 GeV. As such, if the contact 
operator approximation is valid for energies 0(10^ GeV), then the sensitivity of the LHG may be expected to far 
exceed any search of 0(1 — 100) MeV gamma rays. But for such light dark matter, it is easily possible that the 
scale of the mediating particles is very light compared to the energy scale of the LHG. If M* is the mass of a single 
mediating particle exchanged in the s channel, then one roughly finds M ~ M^/y/gxQq, where gx and gq are the 
coupling of the mediator to dark matter and hrst-generation quarks, respectively. As a result, the mediator M* can 
be 0(100 MeV), provided the coupling gq is sufficiently small. In that case, standard “monoanything” LHG searches 
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would be unconstraining, and it is not clear if there will be any signal of new physics at the LHC. 

In particular, if the dark matter annihilation process XX —>■ 771 ° proceeds in the s channel with a cross section 
~ 10“^ pb, then the mediator could be as light as ^ 300 MeV with gx ~ 1, provided gq ^ 10“^. It would be very 
challenging for the LHC to find evidence for such a weakly coupled light mediator above the background of electroweak 
interactions, which can produce missing transverse momentum via neutrinos. 


VIII. CONCLUSIONS 

In this work, we have considered the indirect-detection signals which may be seen by gamma-ray telescopes, sensitive 
to the range 0(1 — 100) MeV for models where light dark matter couples to first-generation quarks. For dark matter 
with mass less than TO,n.± ( 2 m,n.±), dark matter annihilation (decay) is constrained by kinematics to yield very simple 
final states involving pions and photons, resulting in a very distinctive photon signature. 

For the mass range we consider, dark matter s-wave annihilation with {av) ^ 1 pb is already ruled out by constraints 
from Planck, as well as by constraints from COMPTEL, EGRET, and Fermi. But if the process by which the dark 
matter density is generated is nonthermal, then the dark matter annihilation cross section may be much smaller and 
may be probed by future experiments, such as ACT, GRIPS, and GAMMA-400. The most promising strategy is 
observing dwarf spheroidal galaxies, due to the possibility of data-driven background subtraction. A search for diffuse 
emission of photons from dark matter s-wave annihilation may also yield new constraints which are better than those 
currently available, but the sensitivity of this strategy is limited by one’s ability to accurately model the astrophysical 
background. However, future telescopes aiming for better sensitivities are expected to—by construction—more tightly 
constrain dark matter p-wave annihilation, for which current bounds are set by existing diffuse searches. 

If dark matter decays, then one finds, as expected, that dwarf spheroidal galaxies are not as attractive of a target, 
because the decay rate only scales as the density. In this case, bounds from diffuse emission searches may be comparable 
to those from dwarf spheroidal searches, provided one can make some reasonable estimate for the contribution to the 
photon flux arising from astrophysical backgrounds. 

For annihilation (decay) to the 777 ° or 7 r° 7 r° channels, upcoming telescopes could be sensitive to models with 
(av) ^ 10“^ pb (r ^ 10^® s). Sensitivity in the 77 channel could be up to an order of magnitude better. These 
represent vast improvements (up to a few orders of magnitude) over current sensitivity from gamma-ray telescopes 
and from cosmological observations. For dark matter models with mx < 2to^± , current direct-detection experiments 
have limited sensitivity, while the sensitivity of collider searches are highly dependent on the details of the coupling 
of dark matter to the SM. On the other hand, indirect searches for gamma rays arising from dark matter annihilation 
or decay can be ideal probes of this class of models, which generically provide striking photon spectra. This result 
motivates the development of future telescopes that will fill the “MeV gap” in sensitivity and will be able to robustly 
test these models. 
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